Ultracold atomic gas in non-Abelian gauge potentials: the quantum Hall effect supremacy 
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Nowadays it is experimentally feasible to create artificial, and in particular, non-Abelian gauge potentials for 
ultracold atoms trapped in optical lattices. Motivated by this fact, we investigate the fundamental properties of 
an ultracold Fermi gas in a non-Abelian U(2) gauge potential characterized by a constant Wilson loop. Under 
this specific condition, the energy spectrum exhibits a robust band structure with large gaps and reveals a new 
fractal figure. The transverse conductivity is related to topological invariants and is shown to be quantized when 
the Fermi energy lies inside a gap of the spectrum. We demonstrate that the analogue of the integer quantum 
Hall effect for neutral atoms survives the non-Abelian coupling and leads to a striking fractal phase diagram. 
Moreover, this coupling induces an anomalous Hall effect as observed in graphene. 

PACS numbers: 03.75.Lm,67.85.Lm,73.43.-f 



I. INTRODUCTION 

Ultracold atoms in optical lattices offer unprecedented pos- 
sibilities of controlling quantum matter and mimicking the 
systems of condensed-matter and high-energy physics |[TJ O . 
Particularly fascinating is the possibility to study ultracold 
atoms under the influence of strong artificial Abelian and 
non-Abelian "magnetic" fields. The experimental realization 
of artificial Abelian "magnetic" fields, which reproduce the 
physics of electrons in strong magnetic fields, is currently 
achieved through diverse schemes: for atoms in a trap the sim- 
plest way is to rotate the trap lO, while for atoms in optical 
lattices this can be accomplished by combining laser-assisted 
tunneling and lattice acceleration methods |[3l|4l[5l|6l[7l[8l, by 
the means of lattice rotations I^^QMl. or by the immersion 
of atoms in a lattice within a rotating Bose-Einstein conden- 
sate (BEC) fV2\ . Several phenomena were predicted to occur 
in these arrangements such as the Hofstadter "butterfly" 1 13] 
and the "Escher staircase" |4 | in single-particle spectra, vor- 
tex formation ElSllMl, quantum HaU effects |[5lfTQlfT5l[T6l. 
as well as other quantum correlated liquids \17]. 

As shown by one of us in Ref. |18 |, it is simple to gener- 
alize the scheme of Jaksch and Zoller for generating artificial 
Abelian "magnetic" fields fT] in order to mimic artificial non- 
Abelian "magnetic" fields. To this aim we have to consider 
atoms with more internal states ("flavors"). The gauge poten- 
tials that can be realized using standard manipulations, such 
as laser-assisted tunneling and lattice acceleration, can have 
practically arbitrary matrix form in the space of "flavors". In 
such non-Abelian potentials, the single-particle spectrum gen- 
erally depicts a complex structure termed by one of us Hofs- 
tadter "moth" 1 18 1, which is characterized by numerous ex- 
tremely small gaps. The model of Ref. 1 18] has stimulated 
further investigations, including studies of nontrivial quan- 
tum transport properties 1 19 |, as well as studies of the inte- 
ger quantum Hall effect (IQHE) for cold atoms 1 16 |, spatial 
patterns in optical lattices [14J, modifications of the Landau 
levels L20jL and quantum atom optics L21n22ii . 



One should note, however, that the U{2) gauge potentials 
proposed in Ref. [18 1 and used in most of the following works 
are characterized by non-constant Wilson loops: atoms per- 
forming a loop around a plaquette undergo a unitary trans- 
formation which depends linearly on one of the spatial coor- 
dinates. Although such gauge potentials are interesting per 
se, the features characterizing the Hofstadter "moth" result 
from this linear spatial dependence of the Wilson loop, rather 
than from their non-Abelian nature. Indeed, the Hofstadter 
"moth"-like spectrum may actually be found in the standard 
Abelian case with a Wilson loop proportional to x (see Fig. 
1). 

Two of us have shown that cold fermionic atoms trapped 
in optical lattices and subjected to artificial "magnetic" fields 
should exhibit an IQHE |16|. If a static force is applied to 
atoms, for instance by accelerating the lattice, the transverse 
Hall conductivity gives the relation between this external forc- 
ing and the transverse atomic current through the lattice. It 
has been shown that this transverse conductivity is quantized, 
o'xy = — where C is an integer and h is Planck's con- 
stant. Note that this quantity can be easily measured from 
density profiles, as shown recently by Umucalilar et al. |[23ll . 
The quantization of cr^y occurs, however, only if the Fermi en- 
ergy of the system is located inside a gap of the single-particle 
spectrum. While the observation of the IQHE seems to be ex- 
perimentally feasible in Abelian "magnetic" fields, it is hardly 
so in the deeply non-Abelian regime in which the gaps of the 
"moth" become very small 1241 . 

The question therefore arises whether the consideration of 
non-Abelian gauge potentials characterized by a constant Wil- 
son loop could stabilize the spectral gaps and guarantee the ro- 
bustness of the IQHE in ultracold fermionic gases and whether 
an anomalous IQHE, as observed in graphene, can exist in 
such systems. 



In this work we provide affirmative answers to these ques- 
tions by considering the IQHE in a system that features a 
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Figure 1: Energy spectrum E — E(^) in the case of the AbeHan 
gauge potential A = (0, 27r$m^, 0) with x — ma, corresponding 
to the W(9w-C6>«5to? Wilson loop W{m) = g^27r$(i+2m) Compare 
with the Hofstadter "moth" depicted in Fig. 1 of Ref. fl61 or in Fig. 
4ofRef. [25_|. 

non-Abelian gauge potential characterized by specific non- 
commutating constant components and by a constant Wil- 
son loop. We calculate the energy spectrum and we obtain 
a robust band structure with well developed gaps, which dif- 
fers drastically from the case of the gauge potential of Ref. 
ifTSl . In particular, we note the existence of van Hove sin- 
gularities in the density of states and we obtain their analyt- 
ical expression. We then evaluate the conductivity (Jxy for 
neutral currents using topological methods: we express cjxy 
in terms of the topologically invariant Chern numbers associ- 
ated to each energy band 1 26 1 . We eventually present a salient 
fractal phase diagram which represents the integer values of 
the transverse conductivity inside the infinitely many gaps of 
the spectrum. In this way, we show that the IQHE survives in 
the non-Abelian regime, but undergoes strong modifications 
with striking similarity to the anomalous IQHE in graphene 
||27]| : the transverse conductivity suddenly changes sign due 
to the presence of van Hove singularities and is, under certain 
conditions, anomalous because of conical energy spectra. 

II. OPTICAL LATTICE COUPLED TO A NON-ABELIAN 
GAUGE POTENTIAL 

We consider a system of non-interacting two-component 
fermionic atoms trapped in a 2D optical square lattice of unit 
length a, with sites at {x = ma^y = na), with n, m inte- 
gers. The non-interacting limit can be reached using Feshbach 
resonances, or simply at low densities. The optical potential 
is strong, so that the tight-binding approximation holds. The 



Schrodinger equation for a single particle subjected to an arti- 
ficial gauge potential then reads 

V^m— l,n) 

m,n 7 (1) 

where Ux (resp. Uy) is the tunneling operator and ta (resp. 
tb) is the tunneling amplitude in the x (resp. y) direction. In 
the following, we use a as the length, and ta = h = t sls the 
energy units, and set = c = e = 1, except otherwise stated. 
The tunneling operators are related to the gauge potential ac- 
cording to Ux = e*^^ . 

Here we consider a general non-Abelian gauge potential 

A = (aay , 27T^m + (3a x , O) , (2) 

where a and P are parameters, {(Jx^cTy) are Pauli matrices and 
^ is the number of Abelian magnetic flux quanta per unit cell. 

In order to realize such a potential we may consider the 
method of Ref. ifTSl . However, the specific form of this gauge 
potential allows us to consider an even more practical scheme 
based on a generalization of the method currently developed 
by Klein and Jaksch 1121. We may use ^^K atoms in F = 9/2 
or F = 7/2 hyperfine manifolds, or ^Li with F = 1/2. For 
^^K one should optically pump and restrict the atomic dynam- 
ics to the two lowest Zeeman sublevels in each of the hyper- 
fine manifolds. One can then employ different lattice tiltings 
in the x and y directions to perform laser (Raman assisted) 
tunnelings that change the internal states of the atoms; this al- 
lows to control the parameters a and (3, and fixes the tunneling 
rate. Finally, the immersion of the system in a rotating BEG 
will allow to control <l> |12|. In experiments, one routinely 
reaches the values of (laser assisted, or direct) tunneling rates 
in the range of 5-10 kHz (c^ 0.5 /iK), Fermi temperatures of 
the same order, and temperatures T :^ 0.2 Tp ^ 50-100 nK 
(see for instance |2, 3 |). 

The tunneling operators are 2 x 2 unitary matrices, 

Ux = cos a + iay sin a, 

Uy{m)= e^^^*^^ (cos P + iax sin p) , (3) 

which act on the two-component wave function i/jm^n- 

The single-particle Hamiltonian is invariant under trans- 
lations defined by the operators i^m^n = i^m-\-q,n and 
Ty ipm,n = V^m,n+i undcr the condition that <l> = |, where 
p and q are integers. Consequently, the system is restricted to 
a g X 1 super-cell and one can express the wave function as 
V^m,n = e^^'^^e^^y'^Um , with Um a g'-periodic function. The 
wave vector k belongs to the first Brillouin zone, a 2-torus 
defined as /c^; G [0, ^] and ky G [0,27r]. The Schrodinger 
equation ([T]) then reduces to a generalized Harper equation 



3 



^ / cos a sma\ ^i. , /cos a — sma\ 

\^-sma cos ay "^^^ \^sma cos a J ^ 

_^2( cos(27r<l>m + /c^) cos/3 — sin(27r$m + /c^) sin/5\ ^ 
sin(27r<l>m + A:^;) sin/3 cos(27r$m + /c^) cos/3 y 



III. THE NON-ABELIAN REGIME 

Artificial gauge potentials generally induce the following 
non-trivial unitary transformation for atoms hopping around a 
plaquette of the lattice: 

U = UM^^^Mul{m). (5) 



If one sets a = dir or f3 = dir, where d is an integer, the 
matrix U = exp(z27r<l>) is proportional to the identity and the 
system behaves similarly to the Hofstadter model 1 13 1. When 
a = p = {2d^ l)7r/2, where d e Z, one finds that U = 
— exp(z27r<l>) and the system is equivalent to the 7r-flux model 
in which half a flux quanta is added in each plaquette 1271 . In 
these particular cases where U = ± e*^^*^, the system is in the 
Abelian regime. For any other values of the parameters a and 
/3, the matrix /7 is a non-trivial U (2) matrix and the system is 
non- Abelian. 

A gauge invariant quantity which characterizes the system is 
given by the Wilson loop 

W = trU^Uyim + l)U^^Ul{m) 

= 2 e^^^^ (cos^ a + cos 2/3 sin^ a) (7) 

It is straightforward to verify that the system is non- Abelian 
when \W\ 7^ 2, and that W{a,(3) = W{(3,a). In Fig. [2j 
where we show the Wilson loop's magnitude as a function of 
the parameters, \W\ = /3)|, we can easily identify the 

regions corresponding to the Abelian (\W\ = 2) and to the 
non- Abelian regimes (\W\ 7^ 2). We note that the Abelian 
TT-flux regime is reached at a singular point, a = P = {2d + 
l)7r/2, where d G Z. 

We also point out that the statement according to which 
the non- Abelian regime is reached when [Ux^ Uy] ^ 0, and 
which can be found in previous works |[T6l ITSl , is incorrect: 
for the situation where a = /3 = (2(i+l)7r/2, one finds that 
[Uxj Uy] = 2i e^^^^^cFz, while the system is Abelian because 
of its trivial Wilson loop, = | - 2e^'''^\ = 2. 

Contrary to the non- Abelian systems considered in previous 
works (ISIITSKTJI, we emphasize that the gauge potential Eq. 
([2]) leads to a Wilson loop which does not depend on the spa- 
tial coordinates. In the following section, we show that this 
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In the presence of the gauge potential Eq. ([2]), atoms perform- 
ing a loop around a plaquette undergo the unitary transforma- 
tion: 



(6) 
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Figure 2: Wilson loop's magnitude as a function of the parameters 
\W\ — \W{a^(3)\. The system is equivalent to the Abelian Hof- 
stadter model along the lines a = c/tt or /3 = dn, where d G Z 
and is equivalent to the Abelian vr-flux model at singular points, 
a = [3 = (2d+l)7r/2. For any other values of the parameters 
a and (3, the system is non- Abelian. 



feature leads to energy spectra and fractal structures which 
significantly differ from the Hofstadter "moth" lITSl . 



jj _ gi27r<^> ( cos a + cos 2(3 sin + f sin 2a sin 2(3 

— sin2asin^ /3 — isin^ asin2/3 cos^ a 



sin 2a sin^ (3 — i sin^ a sin 2(3 

cos 2/3 sin^ a — | sin 2a sin 2(3 



4 



IV. THE ENERGY SPECTRUM 

The energy spectrum can be obtained through direct diago- 
naHzation ofEq. 

In the AbeHan regime corresponding to a = dir or (3 = dn, 
where d e Z, one finds q doubly-degenerated bands for 
^ = I . In this particular case, the representation of the spec- 
trum as a function of the flux ^ leads to the fractal Hofstadter 
"butterfly" (TSj. For the other Abelian case a = P = |, the 
system behaves according to the 7r-flux lattice: the spectrum 
E = E{^) depicts a Hofstadter "butterfly" which is contained 
between ^ = [0.5; 1.5], i. e. shifted by <l> = 0.5 with re- 
spect to the original "butterfly", and the system remarkably 
describes zero-mass Dirac particles |27|. 

In the non- Abelian regime, which is reached for arbitrary 
values of the parameters (a, the spectrum is constituted of 
2q separated bands as illustrated in Fig. |3] For these general 
situations, the representation of the spectrum as a function of 
the flux $ leads to new interesting features. As in the Abelian 
case, one observes repetitions of similar structures at various 
scales. However, new patterns arise in the non- Abelian case, 
as illustrated in Fig. |4]for a = /3 = J and in Fig. |5]for a = 1 
and (3 = 2. It is worth noticing that for arbitrary values of 
the parameters (a, the spectra show well-developed gaps 
contrasting with the Hofstadter "moth" which appears in the 
non- Abelian system proposed in Ref. |[T8l . We further no- 
tice that the spectrum is periodic with period = 1 and is 
symmetric with respect to £^ = and $ = 0.5. 

In the non- Abelian regime close to = 7r/2, one ob- 
serves that conical intersections are preserved in the energy 
spectrum. As shown in the next section, the particles behave 
similarly to Dirac particles in this non- Abelian region and the 
system exhibits an anomalous quantum Hall effect. 

We eventually note that when the flux ^ = 0, the density 
of states reveals several van Hove singularities at the energies 
E = ±2(1 + cos x) and£^ = ±2(1 — cos x), where x = Ci,p. 
As the flux increases, these singularities evolve and generally 
merge. 



V. INTEGER QUANTUM HALL EFFECT AND THE 
PHASE DIAGRAM 

We evaluate the linear response of the system described 
by Eq. ^ to an external force (lattice acceleration) applied 
along the y direction and we evaluate the transverse conduc- 
tivity using Kubo's formula. Following the method of 
Ref. 1 16], one can generalize the well-known TKNN expres- 
sion Ii28il to the present non- Abelian framework, yielding 

- {dkyUxj\dk^uxj)^dk, (8) 

where uxj is the jth component of the wave function corre- 
sponding to the band Ex such that Hux = Exux, and 




Figure 3: (Color online) Spectrum E — E(kx^ky) for a — (3 — 1 
and $ = I . When the Fermi energy lies within a gap, the transverse 
conductivity h axy is quantized (blue integers). 
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Figure 4: Spectrum E = E{^) for a = 13 = f and $ = gf^, 
where p is an integer. When $ = 0, four van Hove singularities are 
located at ^ = ±(2 + V2) and E = ±{2 - V2). 



refers to the first Brillouin zone of the system. The Fermi en- 
ergy Ep is supposed to lie within a gap of the spectrum. The 
transverse conductivity is then given by the contribution of all 
the states filling the bands Ex < E^ situated below this gap. 

Eq. ^ conceals a profound topological interpretation for 
the transverse conductivity based on the fibre bundle theory 
1291 . In the present framework, such bundles are conceived as 
the product of the parameter space with the non- Abelian 
gauge group U (2). This product space, which is supposed to 
be locally trivial but is generally expected to twist globally, is 
characterized by the non- Abelian Berry's curvature 

^ = {dk^Ay - dk^A^ + [A^,Ay])dk,dky, (9) 

where {A^)ij = {uxi\dk^uxj) is the Berry's connection. The 
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triviality of the fibre bundle is measured by the Chem number 

C{Ex) = / tr^, (10) 

which is a topological invariant and is necessarily an integer. 
Note that each band Ex is associated to a specific fibre bundle, 
on which a Chern number is defined. One eventually finds that 
the Hall-like conductivity Eq. ^ is given by a sum of integer 
Chern numbers, 

As a consequence, the transverse Hall-like conductivity of the 
system evolves by steps corresponding to integer multiples of 
the inverse of Planck's constant and is robust against small 
perturbations. 

The evaluation of these topological invariants leads to a 
complete understanding of the IQHE which takes place in the 
present context. The aim is then to compute the Chern num- 
ber associated to each band Ex of the spectrum. This com- 
putation can be achieved numerically thanks to an efficient 
method developed by Fukui et al [30] and which can be ap- 
plied to our specific system. This method is summarized as 
follows: the Brillouin zone T^, defined by e [0, ^] and 
ky e [0, 27r], is discretized into a lattice constituted by points 
denoted ki = (k^u kyi). On the lattice one defines a curvature 

expressed as 

^i2(/cO = lnUi{ki)U2{ki + i)Ui{ki + 2)-'U2{ki)-\ 

(12) 

where the principal branch of the logarithm with — tt < 
^12 /i < TT is taken, /i is a unit vector in the direction /i, 
and 

U^iki) = Y,{^xjiki)\uxj{ki + fL))/MM^ (13) 

3 

defines a link variable with a normalization factor M^{ki) 
such that \Uij^{ki)\ = 1. The Chern number associated to the 
band Ex is then defined by 

^=iE^-(fc'). (14) 

I 

This method ensures the integral character of the Chern num- 
bers and holds for non-overlapping bands. In the situations 
where the spectrum reveals band crossings, a more general 
definition of the link variables U^{ki) has been proposed in 
Ref. m. 

We first compute the Chern numbers for a specific case, 
illustrated in Fig. [s] For a = f3 = 1 and ^ = |, the 
Chern numbers associated to the six bands are respectively 
1;— 5;6;2;— 5;1. According to Eq. ([TT]), the transverse con- 
ductivity's values associated to the 5 gaps are 1; —4; 2; 4; —1 
as shown in Fig. |3] 

The phase diagram describing the IQHE for our model can 
eventually be drawn. In this diagram we represent the quan- 
tized transverse conductivity as a function of the Fermi en- 
ergy Ef and flux ^. Here we illustrate a representative exam- 
ple of such a phase diagram which was obtained for a = 1, 



(3 = 2 (cf. Fig. [5]). This striking figure differs radically from 
the phase diagrams obtained by Osadchy and Avron in the 
Abelian case ISTl since the Chern numbers associated to the 
gaps are no longer satisfying a simple Diophantine equation 
[32 J . Consequently, the measurement of the transverse con- 
ductivity in this system should show a specific sequence of 
robust plateaus, heralding a new type of quantum Hall effect. 

This new effect is comparable to the IQHE observed in 
Si-MOSFET or the anomalous IQHE observed in graphene 
in the "low flux" regime <l> <C 1 corresponding to experi- 
mentally available magnetic fields. In this regime, the quan- 
tized conductivity evolves monotonically by steps of one be- 
tween sudden changes of sign across the aforementioned van 
Hove singularities (see Fig. |5]). Moreover, in the vicin- 
ity of a, /3 = 7r/2, the quantized conductivity increases by 
double integers because of Dirac points in the energy spec- 
trum, in close similarity with the anomalous IQHE observed 
in graphene. 




-3 -2-1 1 2 3 



Figure 5: (Color online) Spectrum E = E{(^) and phase diagram 
for a = 1,;^ = 2 and ^ = | with q < 97. Warm (resp. cold) 
colors correspond to positive (resp. negative) values of the quantized 
conductivity. Purple corresponds to a null transverse conductivity. 
For $ ^ 1, the quantized conductivity evolves monotonically but 
suddenly changes sign around the van Hove singularities located at 
^ zbl (see the alternation of cold and warm colors). 



VI. CONCLUSIONS 

Summarizing, we have proposed in this paper how to real- 
ize in cold atomic systems a textbook example of non- Abelian 
gauge potential characterized by a constant Wilson loop. Our 
main result is that despite the coupling between the different 
"flavor" components of the single-particle wave functions, the 
spectrum exhibits well-developed gaps of order of 0.1-0.2t, 
i.e. about 50-100 nK. 

The IQHE survives in the deeply non- Abelian regime and 
acquires a unique character specific to the non- Abelian nature 
of the gauge fields. It is characterized by a particular sequence 
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of robust plateaus corresponding to the quantized values of 
the transverse conductivity. Moreover, the non-Abelian cou- 
pling induces controllable van Hove singularities as well as 
an anomalous Hall effect, similar to the effect induced by the 
hexagonal geometry in graphene. Experimental observation 
of this distinctive effect requires to achieve T smaller than the 
gaps, i.e. of order of 10-50 nK, which is demanding but not 
impossible. 

The main experimental challenge consists here in combin- 
ing several established methods into one experiment: laser as- 
sisted tunneling |3 1, BEC immersion 1 12], and density profile 
measurements |23|. 
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